Relevance of phonon dynamics in strongly correlated systems coupled to phonons: 

A Dynamical Mean Field Theory analysis 
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The properties of the electron-phonon interaction in the presence of a sizable electronic repulsion 
at finite doping are studied by investigating the metallic phase of the H[ubbard-Holstein model 
with Dynamical Mean Field Theory. Analyzing the quasiparticle weight at finite doping, we find 
that a large Coulomb repulsion reduces the effect of electron-phonon coupling at low-energy, while 
this reduction is not present at high energy. The renormalization of the electron-phonon coupling 
induced by the Hubbard repulsion depends in a surprisingly strong and non-trivial way on the 
phonon frequency. Our results suggest that phonon might afi^ect differently high-energy and low- 
energy properties and this, together with the efi'ect of phonon dynamics, should be carefully taken 
into account when the effects of the electron-phonon interaction in a strongly correlated system, like 
the superconducting cuprates, are discussed. 
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I. INTRODUCTION 

One of the open problems of high-temperature super- 
conducting cuprates is the role played by the electron- 
phonon interaction. The experimental evidences suggest 
a fairly strong influence of the electron-phonon interac- 
tion on some physical properties, while for some other 
aspects, lattice effects seem to have almost no role. On 
one side, inelastic scattering measurements have shown 
that a specific optical in-plane phonon mode displays 
an anomalously pronounced softeningji^i^ The coupling 
to the same phonon mode has been invoked to explain 
the kink in the nodal electron dispersion detected by 
photoemission,'^-®'^ On the other hand, the almost per- 
fectly linear behavior of the resistivity in a wide range 
of temperatures seems to indicate a little influence of 
phonons on transport properties. Isotope effects also 
show a complex phenomenology: While the the super- 
conducting temperature has a little isotope effect at op- 
timal dopingf^iSiiS the in-plane penetration depth is much 
more sensitive to isotope substitutioniiiiiSii^ Such a large 
(around 5%) effect is usually translated into an isotope 
effect on the effective mass suggesting the presence of 
polaronic carriers in underdoped compounds^ 

The above puzzling scenario leads to a lively de- 
bate which ultimately focuses on whether the effects of 
electron-phonon (e-ph) coupling on different quantities 
are depressed or enhanced by the presence of strong cor- 
relations. Given the intrinsic nonperturbative charac- 
ter of the problem, there is no obvious theoretical ap- 
proach. Different approaches seem indeed to draw con- 
flicting scenarios in studies of the Holstein model for 
the e-ph coupling in the presence of strong correlations, 
described by the Hubbard or the t-J models. Accord- 
ing to Quantum Monte Carlo (QMC) and exact diag- 
onalization (ED) a single hole in the t-J model expe- 



riences an enhanced polaronic effect due to the "pre- 
localizing" mechanism associated to the antiferromag- 
netic spin backgroundii^iiSiiLiS QMC calculations in the 
Hubbard model at finite density and fairly high tem- 
perature suggest that strong correlations favor the small 
transfered momentum electron-phonon vertex (or depress 
it less than its large momentum counterpart), and that 
this quantity increases by increasing repulsion for a win- 
dow of parametersii^ Slave boson approaches suggest 
however that such an effect is just a finite-temperature 
precursor of a phase separation that would take place 
at low temperaturesiSSiSi Such tendency towards phase 
separation, indicated by mean-field approaches also for 
the three-band Hubbard model has been recently con- 
firmed by the more accurate Dynamical Mean Field The- 
ory (DMFT),^ 

Another piece of information conies from the DMFT of 
the half-filled Hubbard-Holstein model in the paramag- 
netic sector (i.e., neglecting the antiferromagnetic order- 
ing), where detailed phase diagrams are available i^^iS^ In 
Ref. |2^we have shown that, close to the Mott transition, 
phonons have different effects on high- and low-energy 
single particle properties (self-energy, spectral weights, 
density of states, . . . ). Namely, high-energy features are 
significantly affected by phonons, whereas the low-energy 
quasiparticle features are basically untouched and they 
coincide with those of an effective purely electronic with 
a slightly weaker repulsion U. In other words, the strong 
correlations strongly reduce the impact of e-ph interac- 
tion on quasiparticle properties, at least when the e-ph 
coupling is not too large. The only residual effect is a 
phonon-induced screening of U which moreover is found 
to vanish linearly with the phonon frequency in the adi- 
abatic limit. 

This variety of results is certainly due to the different 
physical regimes they refer to. In this work we start from 
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our analysis of the half-filled Hubbard-Holstein model 
and relax the half-filling condition, therefore putting our- 
selves away from the Mott transition, in a regime closer 
to that of other approaches. We find that strong corre- 
lations still significantly harm the c-ph interaction. In 
our strongly correlated metal polaronic behavior only es- 
tablishes for e-ph coupling larger by at least a factor two 
than the corresponding values in the absence of strong 
correlations or for a single hole in an antiferromagnetic 
background. Moreover, we find that the interplay of 
electron-electron and e-ph interaction makes the value 
of the phonon frequency quite relevant, particularly for 
its effect on the effective mass. Remarkably, phonon dy- 
namics turns out to be more relevant close to half-filling. 

The paper is organized as follows: in sec. ^ we 
compare the half-filled and finite doping cases of the 
Hubbard-Holstein model. In sec. IIIII we analyze the 
renormalization of the quasiparticlc weight as a func- 
tion of the electron-phonon coupling and the phonon fre- 
quency. In sec. II VI the behavior of the chemical potential 
is studied and the conclusions are then drawn in sec. W\ 



II. DMFT OF THE HUBBARD-HOLSTEIN 
MODEL: HALF-FILLING VS FINITE DOPING 



to Njnax phonon states (ranging from 20 to 40) , and using 
up to iVf, = 9 sites in the conduction bathj^ 

DMFT has been widely use to study the properties of 
the Hubbard model, and a clear framework for the Mott- 
Hubbard transition has been determined. Without en- 
tering the details of these studies, we just recall the qual- 
itative difference between the half-filled system, where 
a metal-insulator occurs at [/ = Uc2 — 3-D, and finite- 
doping systems that have metallic character regardless 
the value of U . Such a distinction makes the effects of 
phonons different in the two cases. 

At half-filling, in the strongly correlated metallic phase 
for U smaller, but not far from Uc2^ the Hubbard model 
presents a clear separation of energy scales, with high- 
energy Hubbard bands well separated from the low- 
energy quasiparticlc peak. In this regime the charge fluc- 
tuations are frozen, therefore strongly harming the e-ph 
interaction, at least with a Holstein coupling. As dis- 
cussed in Ref. the residual effect of phonons can be 
described as a partial screening of the static Hubbard re- 
pulsion, which is ruled by the parameter ljq/U ^ according 
to the expression U^ff = U — rjXD, with rj given by 



The Hamiltonian of the Hubbard-Holstein model reads 



H 



(1) 



where Ci^a (c|g.) and (a|) are, respectively, destruc- 
tion (creation) operators for fermions with spin u and 
for local vibrations of frequency loq on site i, i is the hop- 
ping amplitude, I] is the local Hubbard repulsion and g 
is an electron-phonon coupling constant. In this work we 
will always consider an infinite coordination Bethe-lattice 
with semicircular density of states of semibandwidth D. 
A ~ 2g^ /ujqD is the standard electron-phonon coupling 
and ojqI D is the adiabatic ratio. 

We solve the model by means of Dynamical Mean-Field 
Theory (DMFT), which in recent years has emerged as 
one of the most reliable tools for the analysis of both cor- 
related materials and electron-phonon interactions. The 
method maps the lattice model onto an effective local 
theory which still retains full quantum dynamics, and 
it is therefore expected to be quite accurate for models 
with local interactions as Q). The mean-field correspon- 
dence between the local theory and the original model is 
achieved by imposing a self-consistency condition which 
contains the information about the original lattice. 
In practice, the local problem is described through an 
Anderson-Holstein impurity modelf22iSS whose impurity 
Green's function has to be calculated and used to gen- 
erate a new impurity model through the self-consistency 
equation. To solve the impurity model we use ED,'^'^ 
truncating the infinite phonon Hilbert space allowing up 



2ujo/U 



(2) 



Notice that the above expression for Uef / correctly repro- 
duces the antiadiabatic limit loq/ D — > oo, where the e-ph 
interaction mediates an instantaneous local attraction of 
strength A. The description in terms of an effective Hub- 
bard model with a suitably rescaled repulsion U^f / works 
surprisingly well in the proximity of the Mott transition: 
The full Hubbard-Holstein model and the effective purely 
electronic model display an identical electronic spectrum 
at low energy and, quite remarkably, the same rescal- 
ing for Ueff makes the position of the Hubbard bands 
basically coincide. The effects of the electron-phonon 
coupling close to the Mott transition are well visible in 
the high-energy parts of the electronic spectrum, where 
phonon satellites show up, while the low-energy metallic 
peak turns out to be "protected" by correlation. 

This paper is devoted to the extension of this analysis 
to generic fillings, still limiting ourselves to the paramag- 
netic sector. Conceptually the peculiarity of half-filling is 
the clear separation between high- and low-energy scales, 
which is rapidly lost as the density deviates from one. 
This lack of hierarchy of energy scales makes it more com- 
plicated to draw a simple physical picture such as the one 
described above. On the other hand, the properties of the 
metallic phase at finite doping are not severely affected 
by the antiferromagnetic order, which, if included, would 
instead change considerably the half-filling picture. 



3 



III. QUASIPARTICLE PROPERTIES 

As in Ref. l26l we start our analysis from the inspection 
of the quasiparticle weight 



Z= 1- 



duj 



(3) 



where E' is the real part of the local self-energy E. Due 
to the momentum independence of the self-energy in 
DMFT, Z is also inversely proportional to the quasipar- 
ticle effective mass ratio m* /m. It is therefore clear that 
this quantity measures the metallic nature of the system, 
small values of Z implying poorly metallic situations. 



A. Dependence on electron-phonon coupling 
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At half- filling, in the correlated metal for U < Uc2, 
the quasiparticle weight Z increases with the electron- 
phonon coupling AiS£*2& This counterintuitive behavior, 
in stark contrast with weakly correlated systems, where 
Z decreases with A, is understood in terms of a phonon- 
driven attraction which counteracts the Hubbard repul- 
sion, and it can be quantitatively described within the 
effective picture introduced in sec. ^1^^ By increasing 
A, the effective repulsion decreases, thus the system be- 
comes less correlated, and the value of Z increases. As 
mentioned above, such a reduced effectiveness of the e- 
ph coupling can be associated to the freezing of charge 
fluctuations, to which the phonons are coupled. Close 
to Mott state, most sites are singly occupied, and dou- 
bly occupied and empty sites are a minority. As soon 
as we move away from half- filling, doping, e.g., with 
holes, the number of empty sites proliferates. Therefore 
charge fluctuations are gradually restored, even if they 
are still reduced with respect to a system without elec- 
tron correlations. As a consequence, one can expect a 
stronger electron-phonon signature with respect to half- 
filling. The behavior of Z rapidly becomes more ordi- 
nary, namely Z decreases as a function of A (except, as 
we will see below, for extremely large loq of the order of 
U) meaning that the predominant effect of the coupling 
to the lattice is of localizing nature. 

This is also consistent with a much weaker dependence 
oi Z on U for U ~ 5D (in comparison with the case U ~ 
Uc at half-filling) which would depress the delocalizing 
effect of any phonon-induced variation of U into Uef / . 

In Fig. ^we plot m* /m as a function of A for different 
doping levels and for ujq = 0.2D. In the left panel we 
show the uncorrelated system {U — 0). m*/m is equal 
to 1 for A = 0, it rapidly increases with A, eventually 
reaching a polaronic regime, testified by an exponential 
growth of the effective mass and by the development of 
finite lattice distortions coupled to the electronsi^^i The 
right-hand panel presents instead a strongly correlated 
case(C/ = 5D). Here the value of m* /m at A = 
is strongly dependent on density, since correlations are 



FIG. 1: (Color online) Calculated values of m* /m as a func- 
tion of A for various values of the density. In both cases of 
U — (left) and U = 5D (right) the phonon frequency tJo has 
been fixed to 0.2D. 



more and more effective in localizing the carriers the 
closer we are to half-filling. In practice m* /m ~ 2 already 
for n = 0.50 and it reaches more than 10 for n ~ 0.9. 
The main result is however that the e-ph interaction is 
not able to substantially modify this values determined 
by correlation, up to A = 1.0-1.5 and m*/m is almost fiat 
in this interval, if compared with the left panel. Thus, 
even if doped system do not show the growth of Z with 
A characteristic of half-filling (close to C/c2), the effect of 
e-ph coupling is anyway quantitatively reduced by a siz- 
able amount. A similar information is brought by the 
location of the polaron crossover: If we consider larger 
value of A than in Fig. ^ the e-ph interaction finally 
becomes able to sizably affect the quasiparticle residue 
Z. In the example shown in Fig. [21 (J7 = 5D, n ~ 0.9, 
ujq/D — 0.2), the curve of Z as a function of A clearly dis- 
plays a crossover between a small-A linear behavior and a 
much faster decrease for A > 2-;-2.5. This is precisely the 
signature of a polaron crossover, as shown by many stud- 
ies for the pure Holstein model, where a similar bending 
of the curve occurs for much smaller A ~ 0.8i^2i2^ 

An alternative marker of the polaron crossover 
is the phonon displacement distribution P{X) = 
{^q\X){X\iPo), where j-^o) is the groundstate vector, and 
|Ar)(A'| is the projection operator on the subspace where 
the phonon displacement value X = 1 / \/2MujQ{a + a'') 
(M being the phonon mass) has a given value X. This 
quantity therefore measures the distribution of the local 
distortions. At weak coupling, P{X) has only one peak, 
which only broadens when the coupling increases, but 
when a polaronic groundstate is realized, it presents two 
peaks, corresponding to different distortions associated 
to the different charge statesJ^ The two insets of Fig. 
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FIG. 2; (Color online) Behavior of the quasiparticle weight 
Z (main panel) in the crossover from weak to strong A for 
U — 5D, n = 0.90 and loq — 0.2D. In the two insets the 
phonon distribution function P{X) in the polaron region are 
shown for U = 5D (smaller inset, A — 2.8) and U — (larger 
one, A = 0.875). In both insets n = 0.70 and t^o = 0.2D. The 
light blue shaded area is the full P{X) while the triangles, 
solid circles and dashes represent P{X) restricted to rump = 0, 
1 and 2 respectively. 



El show P{X) for two sets of parameter chosen in order 
to highlight the different way to reahze the bimodality 
in the uncorrelated system (left/larger panel) and in the 
strongly correlated one (right/small panel). The simple 
P{X), shown as a shaded area, does not allow us to dis- 
tinguish the two situations. Therefore we also plot the 
conditioned probability distributions projected onto the 
states in which the impurity has different occupations 
riimp = 0, 1, 2.^^ For C/ = the bimodality is associated 
to a large number of empty and doubly occupied sites 
(that gain more e-ph energy), and a smaller number of 
singly occupied ones. This is a signature of bipolaronic 
groundstate, where pairs of polarons are formed. Since 
we are at finite doping the two distributions have differ- 
ent shapes and heights. For U = 5D the two peaks of 
P{X) are instead associated to empty sites and singly 
occupied ones, since the strong on-site repulsion unfa- 
vors double occupancies. Exploiting the empty sites, the 
system can acquire polaronic groundstate at intermedi- 
ate values of A, while at half- filling this would require to 
completely overcome the Hubbard U. We finally notice 
that in the special half-filling case P{X) is not the ideal 
quantity to identify the polaronic behavior, which, at in- 
termediate values of A, can only appear in the dynamics 
of the excitations (e.g. the behavior of a single hole) 

To summarize the results of this section, we can con- 
clude that for metallic situations with finite doping it 
is not possible to describe the effects of the phonons on 
quasiparticle properties in terms of an effective screened 
Coulomb repulsion, as it happened at half-filling. It is 
anyway still true that the effect of phonons on quasipar- 
ticle properties is substantially weaker than for weakly 
correlated systems, and that the phonon effects can be 



strong only in the high-energy part of the spectra.« The 
polaron crossover is found at values of A which are sen- 
sibly larger than in the absence of strong correlation. 



B. Dependence on phonon frequency 

As it has been discussed in many previous studies, the 
coupling A is not the only parameter which controls the 
properties of e-ph interaction, since also the phonon fre- 
quency plays an important j-ole.'^®''^^"^'* In this section we 
investigate precisely the role of this quantity in strongly 
correlated systems. 

In Fig. 121 we report Z as a function of A, for U — 5D, 
n — 0.70 and different values of ujq. The antiadiabatic 
curve (denoted by loq = oo) is simply obtained for a 
Hubbard model with total repulsion given by [/ — AD. 
In the region A < 1, well before any polaronic behav- 
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FIG. 3: (Color online) Quasiparticle weight Z as a function 
of A for different values of ujo/D. U — 5D and n = 0.70. 

ior, Z is almost linear in A, with a slope strongly de- 
pendent on LjQ. In Fig0| we plot the slope r, defined 
through Z/Zx^Q = 1 + rX. For — >■ oo, where the 
phonons only give rise to an instantaneous attraction, 
which opposes the Hubbard repulsion, Z increases as a 
function of A, and r > 0. This behavior is indeed limited 
to LOo > 5-D, while for smaller frequencies the phonons 
play a more standard role, decreasing Z, as expressed 
by a negative r. More interestingly, r displays a non- 
monotonic behavior by further decreasing ujq. Starting 
from ll!q — 0, it decreases up to loq ~ 0.3Z?, roughly inde- 
pendently on doping, and then rises, remaining negative 
for a wide range of frequencies, and eventually becoming 
positive for quite large ujo/D- While the evolution of r 
from negative to positive can be simply understood in 
terms of a crossover from an adiabatic region, where A 
acts as a localizing force, to an antiadiabatic one where A 
decreases the localizing power of U, the non-monotonic 
behavior suggests that the strongly correlated systems 
displays more relevant energy scales. The same phe- 
nomenon can be observed by plotting the effective mass 
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FIG. 4: (Color online) Coefficient r, defined as the slope of 
Z/Zx^o for small A as a function of uq/D, for n = 0.50, 0.70 
and 0.90. 
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FIG. 5; (Color online) Quasiparticle effective mass vs. ujo for 
A = 0.1, 0.5 and 1.0, at i7 = 5D in a semi-logarithmic plot. 
The three cases n = 0.50, 0.70 and 0.90 are shown on the 
same scale. 



m* /m as a function of ujq at fixed A, as we do in Fig. 
|S1 for the same three values of the doping. The deriva- 
tive of this curve is a measure of the isotope coefficient 
on the effective mass am* ■ For pure e-ph systems m* 
always decreases with luq at fixed A: In the adiabatic 
regime it initially decreases linearly with luq while, for 
ujQ ^ D, it scales as I/ojq^ Therefore the isotope coef- 
ficient Um' = \ ^giogi^o/D) always negative. We find 
instead that, in the presence of strong correlation, m* 
can increase as a function of loq. This effect is larger 
for densities close to the density-driven metal-insulator 
transition: For instance, for n = 0.90, the value of am' 
for ujQ ~ O.ID and A ~ 1.0 is a positive, even if quite 
small, number am* ~ 1 2 x 10~^. For values of ojo 
around 0.3D the effective mass becomes instead almost 
independent on the phonon frequency, then am' changes 
sign and m* /m monotonically approaches the asymptotic 
antiadiabatic value. These results have an interesting 
consequence on the interpretation of the experiments on 



strongly correlated materials: The sign and the magni- 
tude of am* is very sensitive to the values of the electron- 
phonon coupling and of the phonon frequency of the spe- 
cific compoundtSS^" Furthermore, a small value of am' 
can be observed if a compound is in the crossover region 
between the two regimes, but this cannot be interpreted 
as a sign of applicability of the Migdal-Eliashberg theory, 
despite this latter predicts indeed am* — 0. Turning to 
the cuprates, the sizable {am' ^ —0.5) negative isotopic 
effect for the effective mass extracted from the isotopic 
dependence of the penetration depthji^ would suggest 
quite larger values of A than those studied here, putting 
these materials near (or beyond) the polaronic instabil- 
ity. It is well possible that the antiferromagnetic correla- 
tions, neglected in this paper, favor a polaronic behavior 
at least by decreasing doping. Indeed the experimental 
evidences for a polaronic behavior in the cuprates are 
limited to extremely small doping, well inside the anti- 
ferromagnetic phase 

The isotope effect has instead a different behavior at 
half-filhng for U < 1/^2- Close to the Mott-Hubbard 
transition m* decreases with wq, but the absolute value 
of am* is extremely largei^Si^ This can be explained 
once more by the effective- Hubbard picture with Uef / = 
U — rjXD and 77 = 2ujq/U valid at half-filling and for 
U < Uc2- close to the Mott-Hubbard transition, in fact, 
DMFT predicts that, for a pure Hubbard model with 
U = Ueff, m* cx Uc2/{Uc2 - Ueff) i.e. 



U. 



m cx 



c2 



Uc2-U + 2ujoXD/U 



(4) 



An increase in the phonon frequency too at fixed A deter- 
mines then a sizable decrease in the effective mass. The 
final result is oc — 1/Z, in other words the absolute 
value of am* gets larger and larger, the closer one gets 
to the Mott transition. 

Finally we notice that the bare value of r obtained 
in the absence of correlations does not depend strongly 
on the phonon frequency up to wq ~ D. In the adia- 
batic limit by tq — ~N{Ep), where N{Ep) is the non 
interacting density of states at the Fermi energy. For a 
semicircular density of states and n ~ 0.5 1, one gets 
ro ~ —0.58 0.64. Comparing with Fig. 4, it turns out 
that this bare value is larger by a factor 2^6 than the 
corresponding maximum values in the interacting case. 
If we define an effective coupling for the quasiparticles 
according to 1 -I- rA = 1 -I- rgAe//, we see that Ag// is 
rather smaller than A and this reduction is unusually en- 
hanced at small frequency (ro will instead be maximum 
at u!o = 0). At the same time, as discussed in sec. IIIII 
the polaronic crossover is pushed to higher values of A. 



C. Fermi-liquid theory vs DMFT 

One of the main results of the previous analysis is 
that strong correlation tends to reduce the effects of e- 
ph interaction, and the size of this effect depends on the 
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phonon frequency. An important question is whether or 
not the behavior of the quasiparticle weight Z obtained in 
DMFT in the metaUic region away from polaronic insta- 
bihties can be captured, at least quahtatively, by a weak 
couphng approach to the interaction between quasipar- 
ticles and phonons. In particular we can test the valid- 
ity of a Fermi-liquid analysis in which a clear hierarchy 
is supposed: electron-electron correlations create heavy 
quasiparticles, which in turn interact with the phonons 
via a renormalized density vertex. It is understood that 
DMFT is not able to deal with the momentum depen- 
dence of the renormalized vertex. 

We divide the self-energy in two contributions: 



(5) 



where Sa=o is the self-energy of the pure Hubbard model, 
and Sres contains all the additional interaction due to 
bare and correlation-dressed electron-phonon processes. 
We can define a quasiparticle residue Za=o determined 
only by electronic correlations (We recall that E' is the 
real part of E.) 



(6) 



which is related to the full quasiparticle weight Z by 



Z 



1 - ^A=0 



(7) 



Z /Z\=Q may be seen as a wave-function renormalization 
for the quasiparticles created by U due to the additional 
interaction, and the quantity Z^^o^res pl^Y^ the role of 
a phonon-induced self-energy for those quasiparticles. 

In the small-A regime, and assuming that the quasi- 
particle self-energy is linear in w, Eq. |(7J) can be written 
as 



Z 



Z\=a 



1 + rA, 



(8) 



where r coincides with the phenomenological parameter 
introduced in sec. IIII Bl and plotted in Fig. 0] Within 
Fermi-liquid theory we can derive an expression for r at 
the lowest order in A assuming that Y^res only contains 
e-ph processes dressed by U : 



dKe, 

duj 



-N{Ef)*V*A'^ 



-N{Ef)Zx=qA^\D 
(9) 

= N{Ep) /Z\=o is the quasiparticle den- 
V* — XDZ^^Q is the renormalized e-ph 



where NiEp)* 
sity of states , 
interaction between quasiparticles, and A is the vertex 
which couples the electrons to the phonons (a density 
vertex for the Holstein model). Eq. obviously im- 
plies r = —Z\=[)h?N{EF)- This result holds as long as 



one considers only those self-energy diagrams in which 
correlation dresses exclusively the vertex. 

There are only two specific limiting cases in which, 
making use of Ward IdentitieS)22il^ we are able to find 
explicit expressions for the density vertex A: the static 
and the dynamic limit. These read respectively 



and 



ZA=oA(g = 0,w^O) = l 



(10) 



(11) 



where _F'q'^'^-' is the symmetric Landau scattering ampli- 
tude due to electronic processes. Since we are consider- 
ing electronic processes only, we have, in the static limit 
A = K^'^^ / HL^^'^ where k^'^^ is the compressibility of the 
Fermi liquid in the absence of coupling to the lattice, 
while k'q^ is the non-interacting value equal to 2N{Ep). 
Then Eq. ^ can be written, in the static limit, as 



A=0- 



dw 



-N{Ef)Zx=q 



XD 



(12) 



In the opposite dynamic limit, instead, A = 1/Za=o ^^nd 
Eq. becomes 



QUI 



N{Ef) 



z 



XD 



\=o 



(13) 
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FIG. 6: Coefficient r as a function of the density for three 
small values of luq/D. The meaning of the curves labeled as 
static and mixed limit is described in the text. 

We are now in the position to compare the above 
Fermi- liquid analysis with the DMFT results, where all 
the contributions to the electronic self-energy are con- 
sidered. In particular we focus on the small-wo regime, 
where the interplay between attraction and repulsion is 
more subtle, as suggested by the data of Fig. ^ In Fig. 
El we report the coefficient r, for loq/D = 0.2 and for 
two smaller values and compare them, in the range of 
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densities between 0.7 and 0.95 to the static limit result 
r = -N{Ef)Zx=o{h'^''^ /^^'Q'^f of Eq. The dynamic 

limit is not shown because it gives r oc — I/Za^Oj 
thus divergent in the limit n 1. This is evidently not 
comparable with the DMFT results shown in Fig. [SJ 
implying that this limit does not describe the numerical 
results. 

The DMFT results are quite close to the static limit 
(|12|l only for the smallest frequency we considered, but 
they rapidly move away as the frequency rises, and al- 
ready for Wo = 0.2D the difference between the calcula- 
tion and the Fermi-liquid prediction becomes huge. This 
means that a standard Fermi-liquid approach can only 
be applied in the extremely adiabatic regime, and that 
even for small phonon frequencies, there are quantita- 
tively important corrections to the theory. 

As an attempt to provide a reference for the 
intermediate-wo regime, we considered a "mixed limit" 
solution, in which the vertices A are taken as the geomet- 
ric average between the two limiting expressions (|12f) and 
As one can see in Fig. Elthis "mixed limit" solution 
has a behavior similar to the DMFT data. Interestingly, 
this heuristic choice has the advantage of staying finite 
for n ^ 1, while the static limit goes to zero and appears 
to be appropriate only in the extremely small frequency 
limit. 



IV. SCREENING OF COULOMB REPULSION 
BY PHONONS 

Differently from the half-filled case, in the regimes con- 
sidered in this paper we could not estimate Ueff from 
the quasiparticle weight Z since the delocalizing effect of 
the Coulomb screening turns out to be irrelevant. We 
can however discuss how the screening of Coulomb re- 
pulsion due to e-ph interaction influences the behavior of 
the chemical potential /i. 

At half-filling /i is fixed by the particle-hole symmetry 
condition, while now, it needs to be determined in order 
to yield a given value of the density n. 

In Fig. [3 we report /i vs A as obtained within DMFT 
in the Hubbard-Holstein model varying the ratio lO{)/D, 
comparing, for U — 5D, the cases with n — 0.30, 
n = 0.50 and n = 0.90. For all densities fj, decreases 
linearly in A, while the role of lvq changes with the value 
of the density: for very small densities the slope of fi vs 



A monotonically decreases with increasing loq while it in- 
creases for n close to 1 and it is practically constant in ujq 
at quarter filling (ti = 0.50). Such a behavior of the chem- 
ical potential with n can be explained taking a closer look 
to the self-energy I](icj„). The real part of the self-energy 
at zero frequency shifts in fact the chemical potential in 
the interacting case according to ^ = + I]'(zci;„ — s- 0), 
leading to fi\ = fi\=o + S^es(*'^n ^ 0). A first indication 
on the behavior of S(0) and can be derived through 
the Hartree-Fock (HF) approximation. In this light it is 
useful to formally integrate out the phononic degrees of 
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FIG. 7; (Color online) Behavior of the chemical potential as 
a function of A for different values of ujo/D. The cases of 
n = 0.90, n = 0.50 and n = 0.30, all obtained for U = 5D, 
are shown with the same scale. The full and dashed black 
lines are described below in the text. 

freedom in Eq. The result is a retarded phonon- 

mediated interaction 

V{u^) = -XD^^ (14) 

Including the Hubbard repulsion, the total interaction 
between electrons with opposite spin is then given by 
U -\- 2V{u)), while only V{ijj) is present between electrons 
with parallel spin. The Hartree term equals {U — 2XD) ■ 
n/2. The interaction term U + 2V{u!) between electrons 
with spin a and —a contributes in fact with (U —XD)-n/2 
and the remaining —XD ■ n/2 comes from the term V{uj) 
between electrons with the same spin a. 

The equal-spin term gives also rise to a Fock diagram: 



K"''{iLo^) = — / dENiE) -. ^ [f{E) + bi-cuo)] - -. [f{E) + 6(c.o)] , (15) 

2 J \ ILOn — E — UJQ lUJn — E + LOq ' 



where / and b denotes respectively the Fermi and spin. In the small-wg limit this yields the well known 
Bose functions, and N{E) is the density of state per ~XDN{0)iujn result,— and therefore the Hartree-Fock 
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(HF) self-energy at zero frequency is given by the Hartree 
term only, i.e. 

S«^(0) = (;7-2Ai?)f (c^o^O) (16) 

In the opposite limit, i.e. for ^ oo keeping A constant, 
it can be easily seen that Eq. (|15|l is equal to \D ■ [n — 
l)/2, so that one obtains 

E^^(O) = {U^XD) %-^ - oo) (17) 

It is quite obvious that the HF scheme is not ex- 
pected to hold for the large values of the Hubbard re- 
pulsion ([/ = 51?), we are dealing with. For the pure 
Hubbard model it has indeed been found with iterated 
perturbation theory solution of the DMFT, that the 
Hartree- Fock [/-dependence can be canceled by higher 
order contributions. This finding is confirmed by ED- 
DMFT of the pure Hubbard model, where the low fre- 
quency self-energy is almost independent on U for all 
n < lii 

Therefore we can expect that in the antiadiabatic limit, 
where the interaction is exactly a Hubbard one, the whole 
U — XD factor multiplying the density in the HF self- 
energy (see Eq. (|17|) '). is equally counteracted by higher- 
order terms. This is confirmed by our results of Fig. \7\ 
where fi{\) — fJ.\=o ~ ~XD/2 for all the considered den- 
sities in the antiadiabatic regime. This value is nothing 
but the constant term of Eq. (|17|l . which survives the 
cancellation of the density-dependent term. 

In the opposite adiabatic limit, we do not expect the 
phonon part of Hl()|) to follow the same fate of the purely 
electronic one. In this regime, in fact, the phonon- 
mediated attraction is of completely different nature with 
respect to the instantaneous Hubbard repulsion and the 
terms containing A, beyond the HF, are not forced to be- 
have like the Hubbard term. In other words, while for 
large the phonons induce a screening of the Coulomb 
repulsion leading to an effective repulsion Ueff — U—XD, 
at vanishing phonon frequency this effect is not present 
and Ueff — U. On the basis of these considerations, in 
the limit loq/D — s- 0, we expect a cancellation only of 
the [/-term in (fTC|l . leading to /z(A) — ^JL\=o ~ —XDn. 
Such a relation is confirmed by the DMFT data for small 
loq and fits remarkably well with behavior obtained at 
all the considered values of the density, for the chemical 
potential (see the continuous line in Fig. [7|). This also 
clarifies the origin of the opposite behavior displayed by 
jjL for n < 1/2 and n > 1/2. In fact, as it can be seen 
still in Fig. [7| the fully antiadiabatic curve (denoted by 
ujQ = oo)r^ lays below all the other ones for n = 0.30, 
contrary to the case of n = 0.90 in which it lays above. 
Moreover, still in Fig. [7|it can be seen that the chem- 
ical potential is almost completely independent on loq 
for n — 0.50, which is precisely what is predicted by the 
above considerations at quarter filling. This supports our 
idea that the frequency dependence of the A-contribution 
to Id reflects the frequency dependence of the Coulomb 
screening by phonons, even though we cannot put this 



idea on a quantitative basis as we did instead near the 
Mott transition at half-filling. 



V. CONCLUSIONS 

In this paper we have studied the effects of the e-ph 
interaction on a strongly correlated metal for filling dif- 
ferent from n — 1. Choosing a large value of U/D the 
quasiparticle properties of the system are mainly con- 
trolled by the electron-electron correlations and polaronic 
features only appear at values of e-ph coupling A larger 
by at least a factor two than those in the absence of cor- 
relation. However, even if the effective e-ph coupling for 
quasiparticles is reduced by the large Hubbard repulsion, 
we have identified a strong influence of the phonon dy- 
namics. The quasiparticle weight Z, which in DMFT 
is inversely proportional to the effective mass, becomes 
very small because of the localizing effect of the Hub- 
bard repulsion but it is still substantially influenced by 
the value of the phonon frequency and displays an un- 
usual isotope effect. Specifically, since Z depends non 
monotonically on ujq, the isotope coefficient of the effec- 
tive mass changes sign. This change of sign is a strong 
deviation with respect to Migdal-Eliashberg theory, and 
it is not present in the half-filled system, except for pe- 
culiar situations close to the bipolaronic transition and 
for specific values of U/D,^'^^ The effective quasiparti- 
cle e-ph coupling is particularly reduced in the small 
phonon frequency limit. From this point of view our 
DMFT results bear similarities with a mean-field cal- 
culation based on slave bosons and a variational Lang- 
Firsov transformation42r In that approach one obtains 
near half-filling Z/Zx=o ^ 1 - 1.388(a;o/£')A) at large 
U/D and small ujq/D. In the adiabatic limit this ex- 
pression implies a vanishingly small Ae//, while in our 
DMFT Ae// remains finite for ujq ~ 0, as a result of the 
quantum fluctuations that are neglected in the mean- 
field calculation. These values of Ae// are however quite 
smaller than the bare A, in agreement with a Fermi liq- 
uid description in terms of renormalized e-ph vertices in 
the static limit. It is interesting to note that the same 
mean- field approach/citepaolo2 predicts at large U/D a 
critical value Xpoi — Uc/2D for polaron formation which 
is slightly larger than the adiabatic value Xpoi — 1.328 
which is found for spinless fermions in the half-filled Hol- 
stein model and somewhat smaller than our estimates 
Xpoi ~ 2 2.5 at [/ = 5D.^° For a single hole in the t-J 
with J ~ t/3 a much smaller Xpoi — 0.8 is required to 
reach a polaronic behavior. This is understood in terms 
of the pre-localizing mechanism associated with the an- 
tiferromagnetic background. 

Beside the neglect of antiferromagnetic order and spin 
correlations, our DMFT approach is limited by the in- 
ability to capture the /c-dependence of the self-energy, 
and accordingly of the quasiparticle weight. This limi- 
tations does not allow us to describe the momentum de- 
pendence of the effective e-ph coupling that has been 
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proposed to be induced by strong correlations by vari- 
ous authorsii2iS2iSi The momentum dependence can be 
restored by means of cluster extensions of DMFT, such 
as the Cellular Dynamical Mean-Field Theory. It has 
already been shown that, for the purely electronic Hub- 
bard model, a significant momentum dependence estab- 
lishes when the Mott transition is approached already for 
small clusters such as a 2 x 2 plaquettei^ 

A further crucial assumption of the present work is 
the Holstein form of the electron-phonon coupling term. 
Although this kind of description has been justified for 
the cuprates by explicit calculations starting from a full 
three-band modelj^ it must be noted that the symmetry 
of the electronic degrees of freedom which are coupled 
by phonons plays a very important role in the compe- 
tition between the phonon-mediated attraction and the 
Coulomb repulsion. In a model for fuUerenes, in fact, in 
which the electrons are mainly coupled to Jahn- Teller vi- 
brations, it has been found that the attraction between 



electrons is not screened by a strong local repulsion, giv- 
ing rise to unexpected phenomena like high-temperature 
superconductivity driven by strong correlation^^iSSi5& 
More specifically in those models, phonon-mediated su- 
perconductivity is enhanced by the strong correlations 
close to the Mott transition, due to the specific symme- 
try of the phonons, which couple with orbital and spin 
degrees of freedom, as opposed to the Holstein coupling 
with the electron density. 
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